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Correlation between geometry, electronic structure and magnetism of solids is both intriguing and 
elusive. This is particularly strongly manifested in small clusters, where a vast number of unusual 
structures appear. Here, we employ density functional theory in combination with a genetic search 
algorithm, GGA+f/ and a hybrid functional to determine the structure of gas phase Fe^Oy^^ clus¬ 
ters. For Fe^Oy cation clusters we also calculate the corresponding vibration spectra and compare 
them with experiments. We successfully identify FcsOj, Fe40^, Fe40j, Fe^O^ and propose struc¬ 
tures for FceO^. Within the triangular geometric structure of Fe3 0 ^ a non-collinear, ferrimagnetic 
and ferromagnetic state are comparable in energy. Fe40^ and Fe40j are ferrimagnetic with a resid¬ 
ual magnetic moment of 1 due to ionization. Fe^O^ is ferrimagnetic due to the odd number of 
Fe atoms. We compare the electronic structure with bulk magnetite and find Fe40^, Fe40j, FcsO^ 
to be mixed valence clusters. In contrast, in Fe3 0 j and Fe5 0 ^, all Fe are found to be trivalent. 


PACS numbers: 36.40.Cg, 36.40.Mr, 61.46.Bc, 73.22.-f 

In nano technology there is an ever increasing demand 
for increasing the density of electronic and magnetic de¬ 
vices. This continuous downscaling trend drives the in¬ 
terest to electronic and magnetic structures at the atomic 
scale. In essence, two things are required: first, novel 
materials and building blocks with exotic physical prop¬ 
erties. Second, a fundamental knowledge of the physical 
mechanism of magnetism at the sub-nanometer scale. 

Atomic clusters, having highly non-monotonous behav¬ 
ior as a function of size, are a promising model system to 
study the fundamentals of magnetism at the nanoscale 
and below. Such clusters consist of only tens of atoms. 
Quantum mechanics starts to play an essential role at 
this small scale, adding extra degrees of freedom. Since 
these clusters are studied in high vacuum, they are com¬ 
pletely isolated from their environment. 

To use these clusters as a model system, as a starting 
point, a detailed understanding of the relation between 
their geometry and electronic structure is required. 

Even in the bulk, iron oxide has a wide variety of chem¬ 
ical compositions and phases with many interesting phe¬ 
nomena, such as the Verwey transition in magnetite.—^ 

Experiments performed on small gas phase Ee^^O^ clus¬ 
ters beyond the two-atom case are scarce. The structure 
of one and two Ee atoms with oxygen has been studied 
in an argon matrix using infrared spectra.—^ The corre¬ 
sponding vibration frequencies have been identified using 
density functional theory (DET). 

Iron-oxide nanoparticles have been investigated for 
their potential use as catalyst in chemical reactions.— Eur- 
thermore, since the iron-oxygen interaction has a funda¬ 
mental role in many chemical and biological processes, 
there have been quite some studies, both experimen¬ 
tal and theoretical, of the chemical properties of Ee^-O^ 
clusters.—”"— 

The possible coexistence of two structural isomers for 
stoichiometric iron-oxide clusters in the size range n > 5 
was experimentally measured using isomer separation by 
ion mobility mass spectroscopy for Ee^O^ and Ee^O^^i 


{n = 2-9).— Eurthermore, the formation of Ee^-O^ clus¬ 
ters has been studied in the size range {x = 1-52).— 

The number of theoretical studies is, however, man¬ 
ifold. The magic cluster Eei308 was extensively stud¬ 
ied and identified as a cluster with Ci but close to 
point group symmetry.—”"— However, also the geometry 
and electronic structure of other cluster sizes have been 
studied theoretically.—^”"— The prediction of geometric 
structures requires a systematic search of the potential 
energy surface to find the global minimum. 

The majority of theoretical studies were performed 
using Db 1 ^ ^ ^ 1 he number ot works m 

which Ee^O^ clusters were studied with methods beyond 
DET is very limited and restricted to very small cluster 
sizes. Eor EeO^ its reactivity towards H 2 was studied 
on a wave-function-based GASPT2D level.— Eor Ee202 
the molecular and electronic structure were calculated us¬ 
ing both DET and wave-function-based GGSD(T) meth- 
ods and a B 2 U ground state was found.— Eurthermore, 
Ref. [ 25 I reports that B3LYP functional and GGSD(T) 
calculations give the same energy ordering of different 
states, although the energy differences are overestimated 
by the B3LYP approach. 

Recently, the structural evolution of (Ee 203 )^ 
nanoparticles was systematically investigated from the 
Ee 203 cluster towards nano particles with n = 1328.—^ 
In the size range of n = 1-10, an interatomic potential 
was developed and combined with a genetic algorithm in 
search of the lowest-energy isomer. The isomers lowest in 
energy were further optimized using DET and the hybrid 
functional B3LYP. This way, a systematic prediction of 
the cluster structure was done for neutral (Ee 2 03 )^ clus¬ 
ters. 

Because of its high computational burden, in DET the 
geometric structure is often only relaxed into its nearest 
local minimum on the potential energy surface (PES). 
There is no guarantee that this local minimum corre¬ 
sponds to the global minimum. Almost all previous 
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works only consider either random structures or manually 
constructed geometries. However, for increasing cluster 
size these methods become less successful in finding the 
lowest-energy isomer. Genetic algorithms, in which sta¬ 
ble geometries are used to create new structures, proved 
to be efficient in finding the global energy minimum.— 

This method has been successfully used for transition- 
metal oxide clusters.—^ 

Identification of the geometric cluster structure is a 
delicate and computationally demanding task. There¬ 
fore, comparison with an experimental method to confirm 
the theoretical findings is essential. In this work, we com¬ 
bine previously reported experimental vibration spec¬ 
tra:^ with first-principles calculations and a genetic al¬ 
gorithm to determine the geometric structure of cationic 
Fe^^O^ clusters. Of the nine cluster sizes reported in 
Ref. [sO, only the geometric structure of Fe 40 ^ was iden¬ 
tified. In this work, we will also identify the geomet¬ 
ric, electronic, and magnetic structure of Fe 30 ^, Fe 40 ^, 
FesOy and propose structures for FegOj. 


I. COMPUTATIONAL DETAILS 

We employ a genetic algorithm (GA) as is described in 
Ref. [23 in combination with DFT to optimize the cluster 
structures. For this we use the Vienna ah-initio simu¬ 
lation package (VASP)^ using the projector augmented 
wave (PAW) method.—^ Since the geometry optimiza¬ 
tion is the most computationally expensive part of the 
genetic algorithm, we use the PBE+U method^ with 
limited accuracy for the genetic algorithm. For all ob¬ 
tained isomers low in energy, we reoptimized the geo¬ 
metric structure using the hybrid B3LYP functional with 
higher accuracy and consider different magnetic config¬ 
urations. We then calculate the vibration spectra and 
compare them with experimental results. 

Within the DFT framework, functionals based on the 
local density approximation (LDA) or general gradient 
approximation (GGA) fail to describe strongly interact¬ 
ing systems such as transit ion-metal oxides.—^ Due to 
the overestimation of the electron self-interaction, they 
predict metallic behavior instead of the (correct) wide¬ 
band-gap insulator. In an attempt to correct for this 
self-interaction, one can, for example, employ a hybrid 
functional, where a typical amount of 20% of Hartree- 
Fock energy is incorporated into the exchange-correlation 
functional. Especially for the B3LYP functional it has 
been shown that this results in good agreement be¬ 
tween the geometric structure and vibrational spectra 
for clusters . ^ ^ However, hybrid functionals are quite 
computationally expensive compared to LDA and GGA 
functionals. Therefore, in the genetic algorithm we em¬ 
ploy the GGA+U method to take into account that EeO 
clusters are strongly interacting systems. We use the 
rotational invariant implementation introduced by Du- 
darev and a ^lane wave cutoff energy of 300 eV for these 
calculations.— 


The differences between GGA and GGA+U for iron- 
oxide cluster calculations have been analyzed in Ref. Ei 
This study stresses the importance to go beyond GGA for 
transit ion-metal oxide clusters calculations. Aside from 
the well-known difference for the electronic and magnetic 
structure, it even finds a different lowest energy isomer 
than GGA for Ee 32 033 . In our genetic algorithm cal¬ 
culations we use an = 1/ — J of 3 eV for the Ee 
atoms, based on a comparison between B3LYP calcula¬ 
tions and PBE-FU calculations for the smallest cluster, 
Ee304 (see Sec. cm). Eor this comparison we also cal¬ 
culated the mean absolute difference (A) between the 
occupied Kohn-Sham energies (E^) using B3LYP and 
PBE+U: 


A = E 


c^PBE+^y E.B3LYP 

Ei Ej^ 


n 


( 1 ) 


where n is the number of occupied Kohn-Sham levels. 
Note that, the binding distances are only weakly depen¬ 
dent on the used and our value of 3 eV is close to val¬ 
ues used in other works (e.g., 5 eV— , 3.6 eV— , 3.6 eV— ). 

We used the genetic algorithm as described in detail 
in Ref. [13. New geometries are formed by the Deaven- 
Ho cut and splice crossover operation. To determine the 
fitness we used an exponential function. A generation 
typically consists of 20 clusters. It has been shown that 
the geometry of Ee^^O^ clusters only weakly depends on 
the magnetic degree of freedom.— Therefore, we restrict 
ourselves to the ferromagnetic case in our genetic algo¬ 
rithm. 

Eor all obtained isomers low in energy, we reopti¬ 
mized the geometric structure using the hybrid B3LYP 
functional^'^ and consider all possible collinear ori¬ 
entations of the Ee magnetic moments by constraining 
the difference in majority and minority electrons. All 
forces were minimized below eV/A. Standard rec¬ 
ommended PAWs with an energy cutoff of 400.0 eV are 
used. The clusters are placed in a periodic box of a 
size between 11 and 17 A, which we checked to be suf¬ 
ficiently large to eliminate inter cluster interactions for 
each cluster size. Eor the cluster calculations, a single 
/c-point (F) is used. Since we also consider cationic clus¬ 
ters, a positive uniform background charge is added and 
we correct the leading errors in the potential.— All 
simulations were performed without any symmetry con¬ 
straints. The reported symmetry groups are determined 
afterwards within 0.03 A. For the density of states (DOS) 
calculations we used a Gaussian smearing of 0.1 eV for 
visual clarity. 

To obtain the vibration spectra, the Hessian matrix 
of an optimized geometry is calculated by considering 
finite ionic displacements of 0.015 A for all Gartesian co¬ 
ordinates of each atom. The vibration frequencies are 
obtained by diagonalization of the Hessian matrix. The 
absorption intensity is calculated using^*^ 

. 4 , = 974 . 86 ». (^) , ( 2 ) 
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where is the degeneracy of the vibration mode, 
the mass weighted vibrational mode, g the electric dipole 
moment, and 974.86 an empirical factor. A method based 
on four displacements for each ion was also tested but 
yielded the same frequencies and absorption intensities. 
Zero-point vibrational energies (ZPVE) were calculated 
for the isomers lowest in energy of which the vibration 
spectra are also shown. 

For a quantitative comparison between experimen¬ 
tal and calculated vibrational spectra, we calculate the 
Pendry’s reliability factor.— The Pendry’s reliability fac¬ 
tor is a well-established method in low-energy electron 
diffraction (FEED) to quantify the agreement in contin¬ 
uous spectra and has also been applied to vibrational 
spectroscopy.— 

The experimental used infrared multiphoton dissoci¬ 
ation method (IR-MPD) does not only depend on the 
absorption cross section of a vibrational mode, but also 
on the dissociation cross section. Therefore, we use the 
Pendry’s reliability factor to quantify the comparison of 
vibration spectra since it is mainly sensitive to peak posi¬ 
tions opposed to a comparison of squared intensity. This 
peak sensitivity is achieved by comparing the renormal¬ 
ized logarithmic derivative of the intensity 


Y{u) = 


L-\co) 


( 3 ) 


where L{uj) = and W is the typical FWHM of 

the peaks in the spectra. The Pendry’s reliability factor 
is defined as: 



where we integrate over the experimental range of fre¬ 
quencies. Rp values range from 0 to 2, where 0 means 
perfect agreement, 1 uncorrelated spectra, and 2 per¬ 
fect anticorrelation. In practice, Rp values of 0.3 are 
considered acceptable agreement within FEED. Y{uj) is 
strongly dependent on experimental noise and values 
close to zero, hence, we calculate F^xpt(‘^) by fitting the 
experimental spectrum with multiple Forentzian peaks 
and extract the corresponding W. The theoretical fre¬ 
quencies are also convoluted with Forentzian peaks with 
the same W. Rp is always minimized as function of a 
rigid shift of all theoretical frequencies. 

For the calculations on magnetite we used the VASP 
code. We used a Monkhorst grid of 6 x 6 x 2 and an 
energy cutoff of 400 eV. We used the rotationally in¬ 
variant FSDA+I/ implementation by Fichtenstein et al^ 
with effective on-site Coulomb and exchange parameters: 
U = 4.5 eV— and J = 0.89 eV for the Fe ions. 

We used the monoclinic structure as described in 
Refs. [ 3 ^ 14 ^. and calculated the electron density with 56 
atoms in the unit cell. In Ref. [39|, the charge and mag¬ 
netic moment were calculated by integrating the density 
and spin density in a sphere with a radius of 1 A for 


Fe. This radius appears to be chosen such that compara¬ 
ble values with neutron and x-ray diffraction experiments 
were obtained. 

Note, there is no unambiguous way to define these 
radii in systems consisting of two or more atom types. 
Therefore, we checked the correspondence of our results 
to the earlier reported ones and also performed calcula¬ 
tions with a larger radius of 1.3 A for Fe and 0.82 A for O. 
This is a reasonable choice for Fe^O^ clusters since the 
overlap between different spheres is minimal, but most of 
the intra cluster space is covered. 


II. RESUFTS AND DISCUSSION 
A. Magnetite 


Even in the bulk, iron oxide is well known for its wide 
variety of phases and transitions. Magnetite (Fe 304 ), the 
most stable phase of Fe^O^, is for example well known 
for its Verwey transition.—^ Above the transition temper¬ 
ature Tv, the structure is a cubic inverse spinel. Upon 
cooling below Tv, the conductivity decreases by two or¬ 
ders of magnitude due to charge ordering. Furthermore, 
the structure changes to monoclinic. 

Magnetite has the formal chemical formula 
(Fe^^[Fe^^, Fe^^]^ 04 ) where tetrahedral A sites 
are occupied by Fe and B sites contain both divalent 
(Fe^^) and trivalent (Fe^^) iron atoms. Since magnetite 
is a mixed valence system, it is an excellent reference 
system for our cluster calculations to determine their 
valence state and corresponding magnetic moment. 


TABFE I: Spin moments within atomic spheres of 1.3 A for 
the Fe ions in monoclinic Fe 304 . For reference the values 
within a sphere of 1.0 A are also shown. A and B labels are 
consistent with Ref. 1^. 


Site 

Radius sphere 

Spin moment (hb) 
1.3 A 

Spin moment (ps) 

1.0 A 

Fe^+(A) 

-4.02 

-3.78 

Fe^+(B1) 

3.69 

3.45 

Fe®+(B2) 

4.15 

3.93 

Fe®+(B3) 

4.06 

3.84 

Fe^+(B4) 

3.64 

3.40 


In Table HI the spin moments are shown for the dif¬ 
ferent iron ions. The magnetic moments on the A and 
B sites are antiparallel creating a ferrimagnetic struc¬ 
ture. Within the atomic spheres of 1.3 A the Fe^^ and 
Fe^^ ions have a distinct magnetic moment of 4.0 gp 
and 3.7 gp respectively. Note the difference of 0.3 fip is 
much smaller than the 1 jup atomic value and does not 
depend on the size of the atomic sphere used in the range 
between 1.0 and 1.3 A. 
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B. GGA+U 

To determine the optimal in comparison to the 
B3LYP functional for the genetic algorithm, we per¬ 
formed PBE+t/ calculations on the neutral Fe 304 clus¬ 
ter. The results for the electronic DOS are shown in 
Fig. [Hand compared with the hybrid B3LYP functional. 



FIG. 1: (Color online) The density of states for the hybrid 
B3LYP functional and PBF-hF for different values of Feff. 
The average inter atomic distances are shown on the right, 
where Fe-Oi and Fe-02 refer to the Fe-0 distances between 
bridging O atoms (side) and the capping O atom (center), 
respectively. The mean absolute difference A [Eq. [T] between 
the PBE+F and B3LYP energy levels is also shown and is 
minimal for t/eff = 3 eV, indicating the best match in DOS. 

The valence states within -4 and 0 eV are formed by 
hybridized orbitals between the d orbitals of iron and 
the p orbitals of oxygen. For increasing t/, the majority 
spin d orbitals of Fe decrease in energy, whereas HOMO- 
LUMO gap increases. Note that the HOMO-LUMO gap 
of 1.5 eV for = 4 eV still is 0.9 eV smaller than the 
2.4 eV gap for B3LYP. Furthermore, for = 2 and 
3 eV the Fe d DOS features are very similar to those of 
the B3LYP result. To quantify this we also calculated 
the mean absolute difference A [Eq. [1] for the occupied 
levels; the results are shown in Fig. [H A is minimal 
for = 3 eV, indicating the best DOS correspondence 
to B3LYP. We also show the corresponding bonding dis¬ 
tances within the cluster, where Fe-Oi Fe-02 refer 
to the Fe-0 distances between bridging O atoms (side) 
and the capping O atom (center), respectively. Note the 
interatomic distances only change very little with increas¬ 
ing For l/gfj = 3 eV, the binding distances are within 


0.01 A; furthermore, for = 3 eV and B3LYP the oc¬ 
cupied d orbitals of Fe are at comparable energies with re¬ 
spect to the HOMO level. We therefore used = 3 eV 
for our genetic algorithm calculations. 

C. FegO® 

Although the possible number of isomers increases 
rapidly with cluster size, for small systems such as Fe 304 
the number of possibilities is still small. In Fe 304 , the Fe 
atoms can either form a triangle or a chain. For the trian¬ 
gular configuration, two isomers are low in energy. The 
first isomer consists of a ring like structure where the O 
atoms occupy bridging states and one O atom caps the Fe 
triangle as is shown in Fig. Ha). In the second isomer, 
the additional O atom is not located above the center 
but forms an extra bridge between the two ferromagnetic 
(FM) ordered Fe atoms as is shown in Fig. |2](b). 


FeaOS 



FIG. 2: (Golor online) The energy as function of spin mag¬ 
netization for different neutral Fe 304 isomers. The geometric 
figures on the right show the corresponding geometric struc¬ 
ture. O atoms are shown in red, Fe spin up and Fe spin 
down are indicated with orange (red) and green (blue) colors 
(arrows), respectively. For the lowest magnetic states the rel¬ 
ative energy differences are also shown in black. Isomers (a) 
(black line) and (b) (red line) are equally low in energy with 
a ferrimagnetic and ferromagnetic ground state, respectively 
(0 eV). The M = 6 /Xb state of isomer (a) is 14 meV higher 
in energy. 

Figure [2] shows the energy as a function of spin mag¬ 
netic moment for the neutral Fe 304 cluster with four 
different isomers. For all spin magnetizations, the ge¬ 
ometric structure is optimized and shown on the right 
with its magnetic structure lowest in energy. In Fig. [2] 
and the rest of this work, Fe spin up and Fe spin down are 
indicated with orange (red) and green (blue) colors (ar¬ 
rows), respectively. O atoms are shown in red. For the 
neutral cluster, the two triangular isomers are equally 
low in energy with two different magnetic configurations. 
The difference is smaller than 1 meV and therefore be- 
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yond the accuracy of DFT. In isomer (a), as indicated by 
the black line in Fig. [2l the magnetic ground state corre¬ 
sponds to ferromagnetic alignment between the magnetic 
moments on the Fe atoms and a total magnetic moment 
of 14 //b- The Fe-Fe distances are 2.51 A, the Fe-0 dis¬ 
tances for the bridging O atoms and capping O atom are 
1.84 and 1.99 A, respectively. Aside from the FM ground 
state, also the ferrimagnetic state with a spin magneti¬ 
zation of 4 /iB is low in energy and only 14 meV higher 
than the ferromagnetic state. Note we also considered 
a noncollinear magnetic state with M = 0 but this 
magnetic configuration did not turn out to be energeti¬ 
cally stable. 

Isomer (b) is equally low in energy and shown in red 
in Fig. [2j The magnetic ground state corresponds to a 
ferrimagnetic alignment where the two ferromagnetically 
aligned Fe atoms have Fe-O-Fe angles of approximately 
90°. 

We also considered zero point vibrational energies for 
the three lowest-energy levels. When we include these 
into our consideration, the ferromagnetic state, indicated 
by the black line, is lowest in energy, and the M = 4 //b 
and M = 6 /iB states are 17 and 19 meV higher in energy, 
respectively. 


FeaO: 



Spin magnetization [pe] 

FIG. 3: (Color online) Energy of the Fe 30 j isomers as func¬ 
tion of spin magnetization. Figures on the right indicate the 
corresponding structure. The isomer lowest in energy (a) is 
a Fe triangle with three bridge O atoms and one O atom cap¬ 
ping the triangle. For this isomer, the ferrimagnetic 5 /Xb state 
is lowest in energy. The antiferromagnetic 0 /Xb and ferromag¬ 
netic 15 /Xb state are 20 and 58 meV higher in energy, respec¬ 
tively. Note the antiferromagnetic 0 /Xb state corresponds to a 
non-collinear orientation with 120° angles between the spins. 


D. FegOj 

For the cation Fe 304 cluster we also considered ring 
and chain configurations with different oxygen locations. 
For all four isomers we calculated all possible different 
collinear magnetic states. Since an antiferromagnetic 
(AFM) triangle is the most simple example of geometri¬ 
cally frustrated magnetism, we also considered the non- 
collinear state with M = 0 /Xb where all magnetic mo¬ 
ments have 120° angles with respect to each other. The 
results are shown in Fig. [3l For the charged Fe 30 j clus¬ 
ter, the isomer with a Fe triangle where the fourth O 
atom caps the triangle is, like in the neutral cluster, low¬ 
est in energy, as is shown in Fig.HJ Three magnetic states 
are low in energy: 0, 5 and 15 /Xb, with the M = 5 /Xb 
state being lowest in energy, and the non-collinear 0 /Xb 
and ferromagnetic 15 /Xb are 20 meV and 58 meV higher 
in energy respectively. 

The ferrimagnetic state which is lowest in energy, has 
a reduced symmetry (C^) with respect to the ferromag¬ 
netic state and the antiferromagnetic state. This 

could indicate a Jahn-Teller distortion, but could also 
be the result of the inability of DFT to correctly model 
the antiferromagnetic ground state.— ^ However, to dis¬ 
tinguish between these two cases, methods beyond DFT 
such as CASPT2 and CCSD(T) are required and there¬ 
fore beyond the scope of this work. Note that different 
magnetic states only lead to minor differences in the vi¬ 
brational frequencies. 

Interestingly, the typical classical displacement during 
a zero-point vibration in these clusters is of the order 
of 0.03 A. This is of the same order as the typical dif- 



FIG. 4: (Golor online) The neutral (left) and cation (right) 
Fe 304 lowest-energy isomers. Fe spin up and Fe spin down are 
indicated with orange (red) and green (blue) colors (arrows), 
respectively. O atoms are shown in red. The interatomic 
distances are shown in black. The neutral and cation cluster 
have C^y and Cy point group symmetry, respectively. 


ference in inter atomic distances between different mag¬ 
netic states. Therefore, this could lead to interesting phe¬ 
nomena in which, for example, there is a strong coupling 
through exchange between vibrations and magnetism. 

The second triangular isomer of Fe 30 j is 154 meV 
higher in energy and also consists of a ring structure. The 
magnetic state lowest in energy has a magnetic moment 
of 5 /Xb- The Fe-Fe bonding distances are 2.5 and 3.0 A 
between the AFM and FM bonds within the structure. 
The Fe-0 distances vary between 1.7 and 1.9 A. The 
isomer has a C 2 v point group symmetry. 

The third and fourth isomers consist of a linear chain 
of Fe atoms with two O bridging atoms between each Fe 
pair. The two planes can be parallel or perpendicular. 
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FIG. 5: (Color online) The experimental vibration spectra 
of FcsOj and the calculated isomers lowest in energy. The 
reported energy differences include ZPVE. The Pendry’s reli¬ 
ability factor [Eq. 0] is also shown for each isomer. 



EIG. 6: (color online) The neutral (left) and cation (right) 
Ee 405 lowest energy isomers. The neutral cluster has C 2 sym¬ 
metry, whereas the cation cluster has Cs symmetry. 


In Fig. [iKb), we also show the vibration spectrum of 
the ferromagnetic state of this cluster. The Fe-Fe dis¬ 
tances are increased to 2.74 and 3.11 A, respectively. The 
ferromagnetic structure is 514 meV higher in energy. The 
vibration spectrum is similar but slightly shifted to the 
blue due to the increased bonding distances. 


where the latter is lower in energy. Both isomers have a 
magnetic moment of 5 /ip- 

In Fig. [5l both the experimental and calculated vibra¬ 
tion spectra for the different isomers are shown. The 
experimental spectrum consists of three peaks at 540, 
610 and 670 cm~^. The best match is given by isomer 
(a) with calculated vibrations at 505, 630 and 660 cm~^ 
and a corresponding lowest-i?p factor of 0.30, indicat¬ 
ing a reasonable match with the experimental spectrum. 
Since isomer (a) is also the lowest in energy, it is identi¬ 
fied as the experimentally observed structure. 

E. Fe40°^+ 

Fe 405 also consists of a ring structure in which the 
O atoms occupy the bridging sites and one O atom is 
located above the center, as is shown in Fig.[6l The clus¬ 
ter has antiferromagnetic order. However, not all Fe-Fe 
bonds are antiferromagnetic, but also two ferromagnet- 
ically aligned bonds are present. Therefore, the cluster 
has no € 2 ^ point group symmetry but (^ 2 , since Fe-Fe 
and Fe-0 distances vary between 2.72-2.74 A and 1.79- 
2.33 A respectively. The magnetic state with four AFM 
Fe-Fe bonds is 308 meV higher in energy. 

For Fe 40 ^ the isomer lowest in energy consists of the 
same ring structure but is more symmetry broken, since 
the O atom above the ring is off-center as is shown in 
Fig. [6l Therefore the two Fe-Fe distances are 2.69 and 
3.07 A, the Fe-0 distances vary between 1.76 and 2.01 A. 
The isomer has point group symmetry. Two Fe 2 02 
squares are present within the cluster. Isomer (a) has 
a magnetic moment of 1 /ip due to ionization. Interest¬ 
ingly, the ionized cluster has a different magnetic ground 
state with four AFM Fe-Fe bonds opposed to the neutral 
cluster. 





Wavenumber [cm'^j 


FIG. 7: (Golor online) The experimental and calculated vibra¬ 
tion spectra of Fe 40 ^. The isomer shown in (a) is both the 
lowest in energy and Rp [Eq. 0] and can therefore be iden¬ 
tified as the experimentally observed geometrical structure. 
The reported energy differences include ZPVE. 

The second isomer, 459 meV higher in energy, is shown 
in Fig.llKc). This cage-like structure has Cy point group 
symmetry and a magnetic moment of 9 /ip. Figure [7](d) 
shows the third isomer which is 494 meV higher in en¬ 
ergy compared to Fig. [7|(a). The isomer has almost no 
symmetry (C^), and consists of a ring where one Fe-Fe 
bond has two bridging O atoms. The Fe-Fe binding dis¬ 
tances vary between 2.62 and 3.13 A. The isomer has a 
magnetic moment of 1 /ip. 

In the experimental vibration spectrum of Fe 40 ^ 
shown in Fig. [71 five vibration frequencies can be ob- 
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served: 450, 615, 760, 810, and 1070 cm“^. The vibra¬ 
tion at 1070 cni~^ can be identified as a shifted vibration 
in the O 2 messenger attached to the cluster-messenger 
complex and is therefore omitted in the Rp calculation.— 
The best fit is given by isomer Fig. [3(a) with Rp = 0.42, 
which is also the isomer lowest in energy. The calculated 
frequencies: 479, 630, 637, 772 and 796 cm~^ match all 
within 30 cm~^ to the experimental spectrum. Also, the 
relative intensities between different vibrations are very 
similar. Although the ferromagnetic order increases the 
binding distances within the cluster, the changes in the 
vibration spectrum of Fig. [7](b) are small and therefore 
the structure corresponding to Figs. EKa) an d[Zl(b) can be 
identified as the experimentally observed structure and 
the IR-MPD method is not able to resolve the magnetic 
state in this case. 


Fe40r 



Magnetization [|Jb] 


FIG. 9: (color online) Energy as function of magnetization of 
the neutral Fe406 and cationic Fe 40 ^ clusters. The magnetic 
ground state corresponds to a total spin magnetic moment of 
M = 0 and M = 1 /xb for Fe 406 , and Fe 40 j respectively. 


F. Fe40°/+ 

In Ref.—, the Fe 40 ^ cluster was already identified as 
the structure shown in Fig.|8](b). The reported magnetic 
structure was ferrimagnetic with a magnetic moment of 
9 /iB- 






magnetic moments of 5 jip for each Fe atom. Note this 
also matches with an ionic picture in which the Fe atoms 
in Fe 406 have a Fe^^ valence state resulting in an atomic 
magnetic moment of 5 jip. The corresponding structure 
is shown in Fig. [TOl In Ref. [sO is mentioned that the 
symmetry in the M = 10 /is state is reduced from for 
the ferromagnetic state to . In this antiferromagnetic 
ground state, the neutral cluster has D 2 d symmetry. In 
Fe 40 ^ the symmetry is reduced even further to as is 
shown in Fig. [TOl 



FIG. 8: (color online) The experimental and calculated vibra¬ 
tion spectra of Fe 40 ^ for both the previous and new magnetic 
ground state. The vibration frequencies are very similar but 
differ in absorption intensity. The M — 1 state in (a) is 
187 meV lower in energy. 

In our calculations a magnetic state lower in energy 
was found for the same geometric structure for both 
Fe406 and Fe 40 ^. In this state Fe406 and Fe 40 ^ have a 
magnetic moment of 0 and 1 jip respectively as is shown 
in Fig. |9l These structures are 194 and 187 meV lower 
in energy for Fe 406 and Fe 40 ^ in comparison to the 
previously reported state.— The antiferromagnetic mag¬ 
netic ground state of Fe 4 O 0 was also previously reported 
in Ref. [26|. For Fe406 we also calculated a noncollinear 
state where all magnetic moments point towards the cen¬ 
ter of mass, such state with M = 0 is 30 meV higher 
in energy compared to the collinear M = 0 jUp state. 

For the neutral cluster, minima in energy are obtained 
for M = 0, 10, 20 /iB corresponding to flips of atomic 


FIG. 10: (Golor online) The neutral (left) and cation (right) 
Fe406 lowest energy isomers. The neutral cluster has D 2 d 
symmetry, whereas the cation cluster has Cg symmetry. 

Figure IH] shows both calculated and experimental spec¬ 
tra for Fe 40 ^. The vibration spectra for the two calcu¬ 
lated magnetic states in Figs, (a) and [8](b) show very 
similar behavior. The Rp values of isomer Fig. IHKa) 
(0.48) and Fig. IHl(b) (0.39) are both large and indicate a 
better match for isomer Fig.|8](b). Although the spectra 
for Figs. EKa) and Kb) are very similar, the ferrimag¬ 
netic structure has an extra vibration at 720 cm~ with 
small IR absorption. Furthermore, around 550 cm~^, 
vibrations differ slightly in frequency. Since the men¬ 
tioned differences cannot be experimentally resolved, the 
IR-MPD method is unable to resolve between different 
magnetic states and another type of experiments such 
as Stern-Gerlach deflection is required to determine the 
magnetic moment. 
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G. Fe50?/+ 

The neutral Fe 507 cluster has a “basket” geometry 
as is shown in Fig. [TTJ The magnetic ground state is 
ferrimagnetic with a total moment of 4 /ig due to the 
odd number of Fe atoms. The cluster has € 2 ^ symmetry. 



has eight distinct vibrations at 375, 490, 520, 570, 615, 
710, 780, and 830 cm~^ which are best resembled by the 
isomer lowest in energy shown in Fig. Ha), although 
the gap between 615 and 710 cm~^ seems to be underesti¬ 
mated. Note that this also explains the high-Rp factor of 
0.65 for isomer Fig. [T2](a). Similar to Fe 40 ^ and Fe 40 ^ 
the absorption intensities of vibrations in the range of 
300-500 cm~^ are systematically underestimated. The 
individual vibrations of isomer Fig. Ha) are all in agree¬ 
ment within 35 cm~ . Although isomer Fig.[T2](b) has a 
lower Rp = 0.43, the energy difference of 407 meV with 
isomer Fig. Ha) is large and isomer Fig. Hb) has a 
vibration at 450 cm~^ which is not present in the exper¬ 
imental spectrum and lacks the experimental 375 cm~^ 
vibration. Therefore, isomer Fig. [I2](a) can be identified 
as the most probable ground state. 


FIG. 11: (Color online) The neutral (left) and cation (right) 
Fe 507 lowest-energy isomers. The neutral cluster has C 2 V 
symmetry, whereas the cation cluster has no symmetry. 

The cationic structure of Fe 507 is very different and 
shown in Fig. [TTl Like Fe 40 ^, it consists of a cage-like 
structure. The Fe-Fe distances range from 2.7 to 3.1 A. 
Except for the triple bound O atom, all O atoms form 
bridges between two Fe atoms. The ground state has a 
magnetic moment of 5 jUp. The second isomer is similar 


FesO? 



1 ' +1.06 ev¬ 
il , Rp = 0.65 ^ 

- b ' J 

^ 1 1 1^,1 

||j +407 meV] 

W , , 1 

[ a 

1 Rp = 0.65 ] 

II V . . . 1 . . . 11 
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FIG. 12: (Color online) The experimental and calculated vi¬ 
bration spectra of Fe^O^ . The reported energy differences 
include ZPVE. 

to the neutral ’’basket” structure and is 394 meV higher 
in energy as is shown in Eig. Hb). The structure has 
Cg symmetry and a magnetic moment of 5 jUp. However, 
the atomic spin moments have a different arrangement 
for the neutral and cationic state. 

The third isomer is shown in Eig. He) and is 1.04 eV 
higher in energy. It contains two triple bonded O atoms 
and is ferrimagnetic with M = 5 jj^p. 

The experimental vibration spectrum shown in Eig. [12] 


H. FeeOt 

The isomer lowest in energy found for FegOj is shown 
in Eig.[T3land has Cg symmetry where the reflection plane 
is located through Ee atoms 1, 3, and 6. The magnetic 
moment of this isomer is 1 //b- 



FIG. 13: (Color online) The cation FceO^ isomer lowest in 
energy. The cluster has Cg symmetry. 

The second isomer low in energy is shown in Fig. Hb). 
In this isomer no symmetry is present. Compared to the 
lowest found isomer in Fig. Ha) it is 413 meV higher in 
energy and also has a magnetic moment of 1 jUp. 

Figure (HKc) shows the third isomer, which is a dis¬ 
torted octahedral of Fe atoms in which the O atoms cap 
the Fe triangles. The structure is slightly distorted due 
to the AFM order between spins, which lead to slightly 
altered Fe-Fe distances. This isomer is 483 meV higher 
in energy than isomer Fig. Ha). 

Figure [14] also shows the corresponding vibration spec¬ 
tra of the mentioned isomers and the experimental spec¬ 
trum. The experimental spectrum has vibrations at 392, 
420, 500, 730 and 763 cm~^. Note that none of the 
provided isomers match the experimental vibration spec¬ 
trum completely. This is also shown by the large- 
values of 0.56-0.61 for all calculated isomers. The isomer 
lowest in energy Fig. Ha) is the best match since it also 
has vibrations at 420 and 500 cm~ , but the vibrations at 












9 


FeeOj 



Wavenumber [cm'”'] 


FIG. 14: (color online) The experimental and calculated vi¬ 
bration spectra of FceO^. The isomer shown in (a) is the low¬ 
est in energy. The reported energy differences include ZPVE. 


804 and 825 are considerably shifted with respect to 730 
and 763 cni~^. Furthermore, the vibrations at 640, 671, 
and 713 cm~^ are not present in the experimental spec¬ 
trum. The vibration spectra shown in Figs. EKb) and 
He) fit even worse. Therefore, we can not successfully 
identify the FceOj structure. 

Note that our genetic algorithm implementation only 
uses geometry optimization at the DFT level. At cluster 
sizes of Fe 0 O^ and larger, preselection using empirical 
potentials instead of immediate geometry optimization 
using DFT might be more efficient in generating possible 
isomers. 


I. Electronic structure 

In the bulk, iron-oxide materials have many different 
crystal structures such as hematite, wustite, and mag¬ 
netite with all corresponding different electronic struc¬ 
tures. While in hematite only trivalent Fe^^ is present, 
the mixed valence state (Fe^^[Fe^^, Fe^^]^ 04 ) mag¬ 
netite leads to interesting physical phenomena such as 
ferrimagnetic ordering between the sublattices A and 
B and the Verwey transition in which orbital ordering 
leads to a first-order phase transition in the electrical 
conductivity.—^ 

In clusters, stoichiometries corresponding to both 
hematite (Fe 406 ) and magnetite (Fe 304 , FcgOg) and 
other combinations (Fe 405 , Fe 507 ) occur. We therefore 
expect divalent and trivalent Fe cations to be present in 
the reported clusters. There is no unique method to de¬ 
termine the valence state in materials consisting of mul¬ 
tiple types of elements. We therefore compare both the 
local magnetic moments and the local density of states 
(LDOS) for our cluster calculations with bulk magnetite 
results shown in Section III A1 Since the Fe^^ and Fe^^ 
features in the LDOS are very similar for different cluster 


sizes, we show the LDOS of Fe 40 f which contains both 
Fe^^ and Fe^^ in Fig. [151 The LDOS for other cluster 
sizes can be found in the Appendix. 



-8 -6 -4 -2 0 2 4 6 


E-Ehomo [eV] 


FIG. 15: (Golor online) The total, integrated and local den¬ 
sity of states of the Fe atoms for the Fe 40 ^ cluster. The 
trivalent Fe(l), Fe(2) and Fe(3) all show 3d levels at -6 eV 
and small hybridization with O. The divalent Fe(4), however, 
shows strong hybridization and a single level at Ehomo- 


Table El shows the local spin moments of the clus¬ 
ters: Fe 304 , Fe 40 ^, Fe 40 ^, Fe 507 and Fe^Oj. For 
Fe 30 ^ all three Fe atoms have a similar spin moment 
within 0.04 /ig. A comparison with magnetite suggests 
all Fe atoms are trivalent. This agrees with an ionic bond 
model. Furthermore, this is confirmed by the integrated 
and local density of states shown in Appendix |Al The 3d 
peaks around -6 eV correspond to 15 electrons, indicating 
the hybridization between Fe and O is small. Note that, 
the central oxygen atoms 0(4) and 0(7) are partially 
spin polarized. 

For Fe 40 ^, the spin moment of Fe(4) is 0.5 //g lower 
than the other Fe atoms, indicating three trivalent and a 
single divalent atom. The difference is also in agreement 
with the magnetite results. The Fe(4) also breaks the C 2 
symmetry as is shown in Fig. |6l The local (LDOS) and 
integrated density of states are shown in Fig. [151 Note 
that all Fe^^ have 3d peaks around — 6 eV and small 
hybridization with O is present, similar to the Fe 304 
cluster. The LDOS of the divalent Fe(4) atom however 
shows strong hybridization with O and a single minority 
level at Ehomo- 

Whereas Fe406 only contains trivalent Fe,— for Fe 40 ^ 
this is no longer the case due to ionization. As can be 
seen from Table im three trivalent Fe atoms are present, 
together with a single Fe^^ atom. The spin moment is 
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O I 

reduced with respect to Fe , consistent with a higher 
oxidation state than Fe . 

In Fe 507 , only trivalent Fe atoms are present, con¬ 
sistent with an ionic model and the ionized state of the 
cluster. FcgOj, on the other hand, is again a mixed 
valence cluster where the magnetic moment of Fe(4) is 
0.4 yUg lower than the other Fe atoms, indicating Fe(4) 
is divalent. This is also consistent with the LDOS shown 
in Appendix [Al 

Figure [16] shows the density of states for the different 
cationic clusters and magnetite. The calculated band gap 
of 0.2 eV in magnetite is considerably smaller than for 
the reported clusters: around 3 eV for Fe 30 ^ and slightly 
smaller for Fe 40 ^ and Fe 40 ^. Furthermore, whereas 
magnetite has a t 2 g orbital of Fe^^ just below the Fermi 
energy,— in the reported clusters Fe 40 ^ and FcgOj have 
a similar level due to a divalent Fe atom. Note that the 3d 
orbitals of Fe^^ in the clusters are located around 5.5 eV 
below the HOMO level, which is 2 eV higher in energy 
compared to magnetite. 



Fe s 
Fe p 
- Fed 

Op 


E-Ehomo [6V] 


FIG. 16: (Color online) The density of states for Fea,0^ clus¬ 
ters. For these calculations a smearing of 0.15 eV was used 
for convenience of the reader. The HOMO level is located at 
0 eV and the small occupation above the HOMO level is due 
to smearing. 


III. CONCLUSION 


In this work, we have studied the geometric, elec¬ 
tronic and magnetic structure of Fe^^O^ clusters using 
density functional theory. For Fe 304 we compared bind¬ 
ing distances and electronic structure between the hybrid 
B3LYP functional, and different in the PBE+U for¬ 
malism. We found the best match for Ueff = 3 eV. Using 
the PBE+U formalism and a genetic algorithm, many 
possible isomers were considered. Eor isomers low in en¬ 
ergy, all different magnetic configurations were further 
geometrically optimized. Einally, for the cationic clus¬ 
ters we calculated the vibration spectra and compared 
them with experiments to identify the geometric struc¬ 
ture of Ee 30 j, Ee 40 f, Ee 40 ^, Ee 507 and EcgOj. All 
cationic clusters with an even number of Ee atoms have 
a small magnetic moment of 1 /Tb due to ionization. Eur- 
thermore, comparison with bulk magnetite reveals that 
Ee 40 ^, Ee 40 ^ and Ee 0 Oj are mixed valence clusters. 
In contrast, in Ee 304 and Ee 5 07 all Ee are found to be 
trivalent. 
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Appendix A: Local DOS 


In this appendix we show the integrated and local 
DOS of the clusters Fe 304 , Fe 40 ^, Fe 507 , Fe 0 Oj, 
and magnetite. Figures [ITl [181 lEl and [20] show the 
total, integrated and local density of states of Fe 30 j, 
Fe 40 ^, Fe 507 , and Fe 0 Oj, respectively. Of these clus¬ 
ters, Fe 304 and Fe 507 are pure trivalent and the LDOS 
contains 3d peaks at -6 eV and small hybridization be¬ 
tween Fe and O. Fe 40 ^ contains a single tetravalent Fe 
atom, with a similar LDOS compared to Fe^^. The ion¬ 
ized electron is not removed from the 3d levels at -6 eV, 
but from the hybridized levels with oxygen, as can be 
seen from the integrated density of states. Fe 40 j and 
Fe 0 Oj contain a single divalent Fe atom, which has a 
distinct LDOS, in which there are no peaks around -6 eV 
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but strong spin polarized hybridization with oxygen and 
a single occupied minority level at the HOMO level. Even 
in bulk magnetite, as is shown in Fig. [2T1 the same fea¬ 
tures between divalent and trivalent Fe atoms exist. 


FeaO: 
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FIG. 17: (Color online) The total, integrated and local density 
of states of the Fe 30 j cluster. 


FIG. 19: (Color online) The total, integrated, and local den¬ 
sity of states of the Fe^O^ cluster. 



E-Ehomo [6V] 


FIG. 18: (Color online) The total, integrated, and local den¬ 
sity of states of the Fe 40 j cluster. Fe(l) is tetravalent as is 
shown in Tabled 
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TABLE II: The spin moment for Fe^Oy clusters. The atom numbers correspond to the atom numbers shown in Figures [4I6I1QI 
\m and m The spin moment is calculated using atomic spheres of 1.3 and 0.82 A for Fe and O, respectively. 


Cluster 


Spin moment [/i^] 

1 2 

3 

4 

5 

6 

7 

8 

FesOt 

Fe 

-3.84 

3.88 

3.88 







0 

0.56 

0.00 

0.00 

0.22 





Fe 40 ^ 

Fe 

3.89 

-3.84 

3.89 

-3.40 






0 

-0.05 

0.13 

0.20 

-0.05 

0.20 




Fe40et 

Fe 

-3.22 

3.85 

3.85 

-3.79 






0 

0.01 

0.54 

0.01 

-0.25 

0.00 

0.00 



Fe,Ot 

Fe 

3.85 

3.87 

3.89 

-3.83 

-3.80 





0 

0.01 

0.10 

0.03 

0.51 

-0.09 

0.05 

0.12 


Fe^Ot 

Fe 

3.80 

-3.84 

3.85 

-3.47 

-3.84 

3.88 




0 

0.01 

0.51 

0.01 

0.01 

-0.10 

0.17 

-0.10 

0.01 
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FIG. 20: (Color online) The total, integrated, and local den¬ 
sity of states of the FceO^ cluster. All Fe atoms are trivalent 
except for Fe(4), which is divalent. 


Magnetite 



E-Ef [eV] 


FIG. 21: (Color online) The total and local density of states 
of the different Fe atoms in magnetite. The numbering is 
consistent with Table [H Fe^^ and Fe^^ have a similar LDOS 
to clusters although the symmetry is very different. 




















































